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Abstract 



It has been claimed that extreme black holes exhibit a phenomenon of flux 
expulsion for abelian Higgs vortices, irrespective of the relative width of the 
vortex to the black hole. Recent work by two of the authors showed a subtlety 
in the treatment of the event horizon, which cast doubt on this claim. We 
analyse in detail the vortex/extreme black hole system, showing that while 
flux expulsion can occur, it does not do so in all cases. We give analytic proofs 
for both expulsion and penetration of flux, in each case deriving a bound for 
that behaviour. We also present extensive numerical work backing up, and 
refining, these claims, and showing in detail how a vortex can end on a black 
hole in all situations. We also calculate the backreaction of the vortex on the 
geometry, and comment on the more general vortex-black hole system. 
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I. INTRODUCTION 



The story of black hole hair is an ongoing and interesting one. It was thought for some 
time that black holes were relatively bland objects, classified by very few parameters: charge, 
mass, and angular momentum. This picture has changed significantly in the past decade 
with the discovery of various types of solutions which carry other, more exotic, charges — 
such as the colored black holes 1 1 — and solutions with dressed horizons || . What is clear is 
that when a nontrivial topology is allowed for the matter fields, so-called "no-hair" theorems 
can often be evaded. In this paper, we are interested in the question of abelian Higgs hair, 
which may occur because a U(l) vortex can pierce, or even end, on a black hole ||. This 
particular phenomenon is interesting both from the point of view of hair for the black hole 
as well as providing a decay channel for the disintegration of otherwise stable topological 
defects [@-g|. 

Briefly, the results of || showed that is was possible for a vortex solution of the U(l) 
abelian Higgs model — a Nielsen-Olesen vortex || — to thread a Schwarzschild black hole, 
and that the matter fields reacted very little to the presence of the event horizon. Inclusion 
of backreaction of the vortex on the geometry revealed that it was an appropriate smooth 
version of the Aryal-Ford-Vilenkin geometry [|lOj discovered some years previously. Further 
work [0] showed that the conical singularities in other more complicated geometries could 
be smoothed over by the vortex, which allows the exact vacuum instanton for the splitting 
of a conical defect |TT| to be used to construct a smooth instanton for splitting of physical 
topological defects A technical feature of these smooth gravitational instantons was 

that they contained two U(l) fields, the broken U(l) of the abelian Higgs vortex, and 
an unbroken U(l), electromagnetism, needed to give a regular Euclidean section for the 
instanton (although not needed for regularity of the Lorentzian section). Naturally this 
raised the question as to whether the results of [[J were still valid in the presence of this 
extra field. In the papers of Chamblin et. al. [12] it was argued that while the results of 
were qualitatively the same for nonextremal black holes, in the extremal limit a completely 
new phenomenon occurred, and the flux of the vortex was expelled from the black hole, 
rather like flux is expelled from a superconductor. 

The evidence presented in [12] was of the form of analytic arguments for high winding 
vortices, and numerical work representing a black hole/vortex system for a variety of relative 
sizes of black hole to the string in which the flux lines of the vortex appeared to consistently 



wrap the black hole. In a previous comment [13|, two of us pointed out that there were some 



difficulties with the numerical evidence as stated, and that while expulsion was possible for 
thick, or high winding, vortices, it did not appear to occur for thin vortices. This would 
apparently solve a puzzle noted by Chamblin et. al., namely that one could, in principle, 
take a vortex terminating on a near-extreme black hole (see |14| for a detailed discussion of 
selection rules for terminating vortices) and then charge the black hole up to extremality. 
This would appear to be in contradiction with the principle of flux expulsion. However, if 
flux expulsion is not mandatory, then such a puzzle never arises. There is, however, another 
problem if flux is always expelled. It appears that very thin vortices have a higher energy 
wrapping the black hole than if there were no black hole present at all. This would mean 
that the system is unstable and the string expels the black hole, which then raises a physical 
paradox. A string outside a black hole reacts to the black hole's gravitational field, therefore 
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we might expect it to be attracted. On the other hand, the vortex is not charged under the 
massless U(l), therefore it has no reason to feel any repulsion, except for the putative flux 
expulsion force. It is therefore not easy to see what the equilibrium solution of a vortex-black 
hole system would be. Of course this is a rather naive argument, as it ignores any effects 
due to the conical deficit geometry of a gravitating cosmic string, which does produce a 
repulsive force on charges |15[], as well as an attractive force on masses |I6 |, a point which 
we return to in the conclusion after we have explored the issue of gravitational backreaction 
of the vortex. 

In this paper, we submit the abelian Higgs vortex/extreme black hole system to an 
exhaustive analysis, with the intent of pinning down precisely when, or indeed whether, flux 
expulsion can or cannot occur. We give analytic arguments for flux expulsion in certain 
regions of parameter space, and flux penetration in others, and back up the analysis with a 
wealth of numerical data. We discuss the problems of numerical integration of this system 
(and how our work differs from that presented in |T2|]) and how we have avoided these and 
ensured accuracy of the integrations. We also consider in detail the vortex terminating on 
the black hole. We then include a discussion of gravitational backreaction before concluding. 



II. THE ABELIAN HIGGS VORTEX 

We start by reviewing the U(l) vortex in order to establish notation and conventions. 
The action for an abelian Higgs system is 

Si = J d'x^g [D^D^ - \G, V G^ - ±A($t$ - rff] , (1) 

where $ is a complex scalar field, = V M + ieB^ is the usual gauge covariant derivative, 
and Gfj, u the field strength associated with B^. We use units in which h = c = 1 and a 
mostly minus signature. It is conventional to express the field content in a slightly different 
manner in which the physical degrees of freedom are made more manifest by defining real 
fields X, x an d -P M by 

$(x a ) = V X{x a )e ix{xa) (2a) 
B^x a ) = l -[P,(x a ) - V, X (x a )}- (2b) 

These fields represent the physical degrees of freedom of the broken symmetric phase; X is 
the scalar Higgs field, P M the massive vector boson, and x, being a gauge degree of freedom, 
is not a local observable, but can have a globally nontrivial phase factor which indicates 
the presence of a vortex. The existence of vortex solutions in the abelian Higgs model was 
argued by Nielsen and Olesen ||, and in the presence of a vortex § dx = 2ttN, where N is 
the winding number of the vortex. 

In terms of these new variables, the equations of motion are 

Xn 2 

V^WX - P^X + -^-X(X 2 - 1) = (3a) 

y-2 pu 

V M G^ + — - = 0, (3b) 
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where (3 = A/2e 2 is the Bogomol'nyi parameter [17|, and G^ v is the field strength of P h 
The simplest possible vortex solution is that in flat space: 



X = X(R), P fl = NP{R)d fl (f>, (4) 

where R = r^/Xrj, {r, 0} are polar coordinates, and X and P satisfy the coupled second 
order ODE's 

(5a) 
(5b) 

For N = 1, this is the Nielsen-Olesen solution, and gives an isolated vortex for all (3. The 
vortex core consists of two components — a scalar core where the Higgs field differs from 
vacuum, roughly of width 1/y/Xrj, and a gauge core of thickness 0(/3 1//2 /V%7). For higher 
N, the solutions were given in |18[, the principal differences to N = 1 being that the X- 
field is flattened (X ~ R N ) near the core, and the string is correspondingly fattened. An 
additional difference is that for /3 > 1, higher winding strings are unstable to separation into 
N unit winding vortices [|17]]. Figure |l] presents some solutions for (3 — 1. 
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FIG. 1. The Nielsen-Olesen vortex for (3 = 1 and for N = 1 (solid lines), N = 3 (long-dashed 
lines) and N = 10. 
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In this paper however, we are interested in nontrivial solutions in curved space, specif- 
ically in the presence of a charged black hole. This means that our set-up now has three 
length scales, the two string core widths as already mentioned, and the new scale — the black 
hole's outer horizon radius: 

WHiggs ~ m H iggs = ( 6a ) 

w gaug e ~ m^ uge = (6b) 



th = Gm + J G 2 m 2 — Gq 2 (6c) 
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where m is the ADM mass of the black hole, and q its charge (A$ = q/r). We now set 
Xrf = 1, which considerably streamlines our analysis and equations of motion. The gauge 
width of the core has already been replaced by yffi, and we now replace the mass and 
charge of the Reissner-Nordstr0m black hole by M = Gm/wn and Q 2 = Gq 2 /w 2 H . This 
simply means that we have chosen to set the Higgs mass, rather than the Planck mass, to 
unity. There are now no dimensionful quantities, and Newton's constant, G, is a (small) 
number which we will represent by e, where e = 8irGr] 2 . For a GUT string e =O(10~ 6 ), and 
represents the gravitational strength of the string, which will be used in section |VI], when 
we consider the gravitational backreaction of the vortex. 



For now, however, we ignore the gravitational backreaction of the vortex and, as in |]12 
treat the vortex in the background Reissner-Nordstr0m geometry: 



ds 2 



— r 



d9 2 + sin 2 



2M + <a~V 



(7) 



For the moment, it is irrelevant whether the geometry is the result of an electric or magnetic 
potential since we only require the equations for the vortex fields in the Reissner-Nordstr0m 
geometry. Substituting (0) in @, and assuming the form (f|) for the P^-field gives 



2Mr + Q 2 )d r X 



1 XN 2 P 2 
o ■ n do (sin 6d e X) + , 2 
r z sin r A sin V 



\X(X 2 - I) 



1 - 




sin 9 



d e P 
sin 6 



X 2 P 

P 



(8a) 
0. (8b) 



In general, these equations are intractable analytically, however, as in || we can extract 
some information in a particular limit — the 'thin string limit' where we assume M»l; the 
value of Q is irrelevant. First note that if we write R = rsin^, and make the assumption 
that X = X(R), P = P(R), then (§) becomes 



„ P' X 2 P 
-P + — + 

R (3 



R?_ 

2 



R 2 (2M Q 2 



r 



Q 



P" 



' R 



2M 



P" - 



P^ 
R 



(9a) 
(9b) 



However, note that the RHS of each of these equations has the form R 2 /r 2 times terms of 
order unity. Near the core R is order unity, hence the RHS is 0(M -2 ), thus we have the flat 
space equations (f|) for which the solutions are well known, and have the property that away 
from the core the vortex fields tend to their vacuum values exponentially rapidly. Therefore, 
by the time the premultiplying term in the RHS of @ becomes significant, the fields are 
essentially in vacuum, so any corrections will be negligible. We can therefore regard the flat 
space solutions (as functions of R — r sin 6) as a good approximation to the true solution 
in the thin string limit. Note that this form of the solution pierces the horizon and does 
not depend on the value of Q, therefore, using this argument, one would expect that thin 
strings always penetrate the event horizon of a black hole, whether or not it is extremal. 

The argument developed so far starts from the thin string limit, but there is another 
limit in which considerable information can be extracted analytically, and which leads to 
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the expectation that vortices are expelled from extremal horizons. This regime, which can 
be regarded as a 'thick string limit' complementary to the one above, is attained for large 
winding number N. As shown in |18| , when N is large the size of the vortex grows like \^N 



> 

N 



and the unbroken phase inside the core is approached increasingly faster, X ~ R . Consider 
then a black hole that sits well inside the vortex core. There, the field is expected to be very 
close to the symmetric phase, so it seems reasonable to neglect the last term in (|8|b). Then 
the equation can be solved by 

P «1 -p(r 2 -Q 2 )sin 2 6, (10) 

with p « 1/ (2Ny/]3) [|l^]. From here we see that, in the extremal limit in which the horizon 
is at r = Q, the magnetic flux across the horizon, given by Gq^ = dgP, vanishes. Moreover, 
it is possible to solve for the Higgs field X by setting X = [b(r) sin 9} N , and keeping only 
leading terms in 1/N. One finds, near the horizon, 

d(\ogb) oc . = = = =, (11) 

1 S ; y/r* - 2Mr + Q l M T _ r+ )( r _ r _) 1 ; 



so that, if the black hole is extremal (r + = r_) then b ~ r — r+. Hence, X vanishes as well 
on the extremal horizon. Furthermore, a study of the energy of the configuration shows that 



it is favourable for the black hole to remain inside the vortex core |12| . 

Indeed, the behavior of magnetic fields in the vicinity of extremal horizons has been 
studied in more generality in fl9| and the expulsion of the flux — a phenomenon remarkably 



analogous to the Meissner effect in superconductors — has been argued to be generic. Ex- 
tremal horizons tend to repel magnetic fields, at least if the latter are in, or approach, a 
phase of unbroken symmetry, like in the core of the vortex. The argument above for large 
N vortices, however, is not fully conclusive, since it remains the possibility that corrections 
of higher order in M or in 1/N spoil the exact expulsion phenomenon. 



III. ANALYTIC ARGUMENTS AND BOUNDS 



In the previous section we summarized the arguments of [|12j in favour of flux expulsion, 
and extended the thin string arguments of @j which seemed to indicate flux penetration. 
Both of these arguments appear compelling, and we must examine the system closely to see 
what definite information can be extracted. First, note that all methods are in agreement 
that there is flux penetration for nonextremal black holes, therefore for the rest of this 
section, we will only be considering extremal black holes, for which the metric is 



ds 2 



1 



— dt 2 
r J 



1 



—) dr 2 -r 2 d6 2 -r 2 sin 2 Odcf. 
r J 



This gives the vortex field equations 

1 



V - M) 2 d r X 



r 2 sin 6 



v t\t2 p2 

d e (sinedeX) + , + \X(X 2 
r z sm 



r 



^)V 

r / 



sin 9 



d e P 
sin 6 



X 2 P 



o 

0. 



(12) 

(13a) 
(13b) 
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which as before are not analytically soluble, however, we can extract quite a bit of informa- 
tion about the solutions due to the nature of the geometry near r = M. 

The first observation is that if we set r = M, then the equations for the horizon actually 
decouple from the exterior geometry: 



sin 9 



v j\t2 p2 

d e (smOdgX) + . 2Q + \M 2 X[X 2 



sin 2 9 
sin 9dg 



2 

dgP_ 

sin 9 



M 2 X 2 P 







0. 



(14a) 
(14b) 



a phenomenon which does not occur in the nonextremal case. This means that the vortex 
equations on the horizon are now ODE's and therefore easier to handle. Note that the flux 
expulsion solution (X = 0, P = 1) always solves (|14|) , therefore we cannot use any analysis of 
these equations to demonstrate flux expulsion, but we can potentially show the nonexistence 
of a penetration solution. 

Therefore, assume that a piercing solution to the vortex equations does exist throughout 
the spacetime, this means that a piercing solution must exist on the horizon. This requires a 
nontrivial solution (X(9), P(9)) which is symmetric around 9 = n/2 at which point X has a 
maximum and P a minimum. Let X m and P m be the extremal values of X and P attained. 
In addition, expanding (|TJ]) near the poles indicates that Pg = at the poles. Therefore, 
there exists a #0 (< ir/2) such that P^g = at 9 , and P t g(9 ) < 0. 

The structure of the proof is as follows; we use the properties of the solution at n/2 to 
derive an upper bound on P and P t gg there. Then we use the behaviour of P t g to derive a 
lower bound on P <ge . At the very least, these must be consistent for a piercing solution to 
exist, therefore, if the inequalities are incompatible we conclude that the core solution is the 
only solution on the horizon and therefore flux expulsion must occur. 

To provide the upper bound, consider the X equation at 9 = n/2. Since X^g$ < 0, 
equation (|14a|) implies P^ < §4^(1 - X^ 



< \ and hence 



M 2 M 3 



2 \l/2 



< 



y/2 M 3 

373 W 



(15) 



where the final inequality is obtained by maximizing over X m . 

For the lower bound, on the other hand note that at 9 \P$\ takes its largest value, 
\P,e \ = jfX 2 Ptan9 < ^tan# , and hence 

— 1= > Pm > 1 — — -Pfln ^ \P (} 

Assuming M < \/2A, this gives | - 9 < cot 9 < 




tM 2 



2/3 1 



one can show that P^g has a maximum at 7r/2, hence 



V2N 



. But then for M 2 X 2 m < 2(3 



P,g(7l/2)-P fi (9 ) 



7T 



/2-0 o 



> 



nj M 2 



M V 



v 1 V2N) 



(17) 
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Therefore, by comparing (15) and (|17j), we see that an absolute minimum requirement for a 
piercing solution is the consistency of these two bounds, i.e. 



V2n 2 M 5 > 12V3P 2 N (l - . 

Turning this around therefore, and writing Ai = 77^, we may conclude that the vortex 
flux lines must be expelled from an extreme Reissner-Nordstr0m black hole if 

m> < *&g_^£__ (18) 



(1-.M) 2 2vr 2 iV 4 4iV 4 ' 

For N = (3 = 1, this gives M < 0.7; note that this is a rather weak bound, in fact we would 
expect flux expulsion to be mandatory for M somewhat in excess of 0.7, but this method at 
least provides an analytic proof giving a definite bound for M. Numerical work (section fV\ j 
actually places this bound at about 1.9 (see figure |]). 

It is interesting to note the variation of this bound with iV and (3. For large N, Ai 5 < 
0(N~ i ), or M < 0(N 1 ^ 5 ). This means that the larger N, the larger the black hole can be 
and still have flux expulsion. This is in agreement with the observation that large N vortices 
are thicker than their single winding number counterparts, therefore we would expect flux 



expulsion to occur more readily. Indeed, this is how |L2[ originally argued for flux expulsion. 
For (3 — » we see that M < 0(/3 2//5 ), i.e. that flux expulsion only occurs for extremely 

small black holes. To understand this, recall that the fall-off of the P field is P ~ e~ R l^ 
therefore as (3 — > 0, the magnetic flux core of the string is getting smaller as This is 
consistent with the above bound. It is therefore interesting to look at large (3, since in this 
limit, the magnetic core becomes very diffuse, and we are left with the Higgs core, which 
is no longer exponentially cut-off, but follows a power law cut-off approaching that of the 
global string. For large (3, the bound (0) gives M < l-0(/3~ 2 ), hence M < \/2N+0(f3~ 2 ). 
Therefore, for small charge Higgs, or global strings, we expect flux expulsion to occur for M 
of order the winding number of the string. 

Having shown that flux expulsion must occur for sufficiently thick strings, what of the 
argument of the previous section, which appeared to indicate that a thin string would pierce 
the event horizon? Obviously, since the core is always a solution to the horizon system 
of equations, we cannot use an argument based on this system to argue flux penetration, 
but instead we must look at the full PDE system of equations in the exterior region of the 
horizon, equations fll3l). 

Now, assume that there is a flux expulsion solution, then on r = M, X = and P — 1. 
Therefore, near r = M, M 2 X 2 < 1 and [(r - M) 2 X r ] jr > 0; hence 

\M 2 X sin 2 6 + sin 6d e sin 6d e X < XN 2 P 2 < XN 2 (19) 

in this region. Now, we know that X is symmetric around tc/2, peaking at some maximum 
X m , and also that sin^c^X vanishes at 0, tc/2 and 7r; therefore let #0 be the value at which 
dgsmOdgX = 0, which must satisfy |M 2 sin 2 6 , < N 2 . If M < V2N, then this inequality is 
clearly satisfied, so we now take M > V^N, and let a > 9 be defined by M 2 sin 2 a = 2N 2 . 
Integrating (|19|) on the range (#,7r/2), for 9 > a gives 
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X fi {B) >X(6) [|M 2 cot# + A 2 cscfllntanfl/2] . (20) 

But since X$q < on [9o,ir/2], we can deduce X fi [G) < X{d) 9 Zff o) < fg, hence for 
consistency with (EDI) 



^ ><«-<») 



esc 2 a cot 6> + csc6»lntan6»/2 (21) 



must hold over the range 9 e (a, vr/2) for the expulsion solution to hold. The actual bound 
on M is then obtained by plotting these curves and determining for which M this inequality 
is always satisfied. For N = (3 = 1 we find that for M 2 > 8.5, this inequality is violated, 
hence the vortex must pierce the horizon in this case. 

For larger N, the lower bound on M for a piercing solution to be forced does increase, 
however, the ratio M/N actually decreases. The LHS of (|21| ) is N~ 2 and it is easy to see that 
this requires a = ir/2 - 0(N~ 1/2 ). From this, we therefore obtain M > V2N [1 + 0(iV _1 )]. 
Note that this argument does not depend on j3. This means that for large (3 we still get 
piercing solutions for the same range of M. Since f3 —>■ oo corresponds to the global string, 
this should not be too surprising. For f3 — > 0, although we expect some drop in the value of 
M, this method is unable to detect this. 

When the vortex radius is much larger than the black hole radius one can find approx- 
imate explicit solutions for the fields near the horizon. Indeed, one can construct these 
solutions for arbitrary winding number N, therefore generalizing the solutions in jn| . 

The solution for the magnetic field is readily found by noticing that if M (or better, 
Mj sj~f$) is very small, then inside the vortex and close to the horizon the gauge field is well 



approximated by the solution to the massless field equation fli"2|,|T9" 



P ~ 1 - 2Mp(r - M) sin 2 9, (22) 

where p is an integration constant equal to twice the magnetic field strength at the center of 
the core. Next, by looking at the Higgs field equation, we can see that close to the horizon 
the potential term |X(X 2 — 1) is suppressed by a factor M 2 which we are taking to be 
small. Therefore it can be neglected. After setting P w 1, the resulting equation can be 
solved with X a function exclusively of (r — M) sin 9, 

X ~ k [(r-M) sin 6] N , (23) 

where k is another integration constant. We see that on the horizon both the gauge and the 
Higgs field are expelled. Comparison between these approximate solutions and the results 
of the numerical calculations in the next section shows good agreement. 



IV. NUMERICAL RESULTS 



To solve equations @ numerically, we have used the technique developed in 0, which 
consists of relaxing initial configurations of the fields X and P on the (rectangularly) dis- 
cretized plane, (r, 9) — > (r^ = r + + idr, 9j = jd9). We therefore replace the fields by their 
values on this grid, X(r, 9) — > Xij = X(ri,9j) (and similarly for P), and the differential 
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operators by suitably discretized versions. Adopting the notation of M and 



Xoo 

-^00 



X{ y j, X± Q 



00 ) ^ 



00 



-> P n ew , where 



Xij±i), we find that the discretized version of 



(that is, 
1) ^ 



ynew 
A 00 



nnew 

M)o 



2 Am 

r \ r 



-X+o— -Jf- 
2Ar 



+ 



cot 9 Xq+—Xq- 
r 2 2A9 



Xq++Xq- 
Ar 2 



+ 



r 2 A6» 2 



2i\i 



Ar 



( NPon V 
\ r sin 9 J 





- + 




_ 2M _|_ Q__^ 

r r 2 j 


P0++P0- 
Ar 2 


, P 0+ +P - 
1 r 2 Ad 2 




2 
Ar 2 


+ 


o X 2 
r 2 A<9 2 ^ /3 







(24a) 
(24b) 



There is, however, a subtlety in this process: relaxation methods usually require that the 
values of the fields be fixed at all the boundaries of the domain of integration, and although 
we know the asymptotic values of X and P at r — > 00 (the vacuum) and at 9 — > 0, it (the 
string core values), the configuration at the horizon r = r + is in fact the main result we 
expect from this numerical calculation. The solution to this problem conceived in |3j was 
to update the values of the fields at the horizon immediately after updating the interior of 
the grid. Note that this still requires an initial guess for the fields on the horizon — a crucial 
point we will return to later. Replacing r = r + in (0), we obtain equations on the horizon: 



r + — r_ dX 
dr 



r + — r 

r 2 



(sm9d e X) + \x (x 2 - 1 
ri sin 9 2 v 



N 2 XP 2 
rl sin 2 6 



dP 
dr 



sin 9 



do 



r=rj r 



'dep; 

sin 9 



+ 



X 2 P 







(25a) 
(25b) 



clearly, these equations reduce to ( |I~4"D in the extremal case. We discretize this in the same 
way that we discretized the equations on the interior of the grid (except that we must now 
take discretized differential operators that do not depend on X_ or P_o). The resulting 
equations are 



X, 



00 



00 



A 00 



nnew 
M)0 



X0+—X0- 
4A6» 



Ar 



+ Ae 2 + T t (^"00 ■ ' 2 v mi, 9 j 



P0++P0- 

2A0 2 



L > T 2 V sin6> J 

cot^- Po+ - p °- 



4A0 



y/Af 2 -Q 2 
Ar 



(26a) 
(26b) 



A<9 2 + 2/3^00 



The process of updating the interior of the grid and then the horizon at each iteration was 
carried on until the modulus of the largest relative correction on the grid became smaller 
than some e: 



max 







pncw pold 

i,j i,j 


, max 

i,3 


Void 


pold 

i,j 



< e. 



(27) 



(i and j run over the entire grid, including the horizon.) 

The results obtained by our implementation of this method were compared with the 
plots of ||, and we found a satisfactory agreement; for instance, figure ^ shows the contours 
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FIG. 2. Contours of X and P for M = 10, Q = 0, (3 = 1/2, JV 
= Nq = r m = 100. (The dashed semi-ellipse represents the horizon.) 



100, £ 



10" 4 and 



of X and P for M = 10,(3 = 1/2 and N = 100; it can be directly compared with figure 3 
from ||. (r m is the maximum value for r, approximating r — > oo.) 

As in ||, we have found that the thin string approximation is excellent for thin vortices, 
and is a reasonable approximation even for thicker vortices, which tend to pinch slightly 
near the horizon. 

We then turned our attention towards charged black holes, comparing now our results 
with those of Chamblin et. al. [12]. We found, as they did, that in nonextremal cases the 
picture remains qualitatively the same as for uncharged black holes. For extremal black 
holes, however, our results differ from their original claims. As reported in [13|, we find 
that the claimed expulsion of the matter fields for thin strings in this limit is the result of 
a loophole in the numerical method (when applied to extremal black holes), which does not 
take into consideration the decoupling of the horizon from the main grid. [This was shown 
in Eq. (|14|) and can also be seen from (|26|): if Q = M, these equations do not contain the 
terms X +0 or P + o] In this case, the core configuration is always an exact solution on the 
horizon; bearing in mind that the relaxation method updates the fields from a user-supplied 
initial guess, we see that if the initial guess made on the horizon is core, then the horizon 
will never be updated! In fact, [|EJ always started from this guess on the boundary, and 
therefore always obtained wrapping solutions in the extremal case. 

Obviously, this choice of initial conditions is important, moreover, since we are dealing 
with a nonlinear system of PDE's, there is no reason for different initial configurations to 
relax to the same solution. For this reason, we have considered the following three initial 
data sets on the horizon: 

• Core: X = 0,P = 1, 

• Vacuum: X = 1,P = 0, 

• Sinusoidal: X(9) = sin0,P= 1. 

(The sinusoidal guess was chosen because it interpolates smoothly, and in a simple and 
convenient way, the strings attached to the North and South poles of the black hole.) 

Figure [5] displays the solutions obtained from the three initial configurations. By com- 
puting and comparing the energy densities and total energies of the fields on the grid, we 
were able to determine that inside the grid the three solutions were identical; on the horizon, 
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however, the solution relaxed from the core guess was a string core, whereas the solutions 
obtained from the vacuum and the sinusoidal guesses were both the vacuum. The compari- 
son also showed that the wrapping solution had a higher total energy than the piercing one; 
the difference is of course due to the jump of the fields from the horizon to the interior of 
the grid. (Making the stepsize dr smaller failed to smooth out this sharp jump.) 

To summarize, we have three (physical and numerical) reasons to prefer the piercing 
solution to the wrapping one for thin strings in the extremal limit: it is smooth, numerically 
more robust, and energetically favourable. 
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FIG. 3. Contours of X and P for a core guess (top) and a vacuum (or sinusoidal) initial guess. 
The parameters are M = Q = 10, j3 = N = 1, e = 1CT 4 , N r = N e = 100. 



To determine how the transition from a piercing to a wrapping solution occurs as we 
thicken the string, we take advantage of the fact that, on the horizon, we now have ODE's. 
This allows for much quicker and more accurate numerical methods; for the following cal- 
culations we have used the relaxation routine solvde of Ref. |2(| chapter 17], starting from 
a vacuum guess. 

Figure |] shows the solutions on the horizon for /3 — N — 1. For massive black holes 
(or, equivalently, thin strings), the fields adopt a vacuum profile on most of the horizon 
(symmetrically about the equator 6 = tt/2) and interpolate smoothly to their fixed core 
values on the poles. The shapes of X and P remain the same for all values of M: X 
has non-vanishing ^-derivatives at the poles, and a maximum at the equator; P has zero 
derivatives at the poles, and a minimum at the equator. As we thicken the string, the 
maximum of X and the minimum of P move away from the vacuum values. For N = (3 = 1 
this transition is gentle at first, but accelerates suddenly, as if the string had crossed a critical 
width at which it is not able to pierce the horizon any more (see figure ||] for an example of 
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expelled solution). This can be observed on figure |6], which shows the evolution of X(jr/2) 
and P(7r/2) with the inverse energy of the black hole for a variety of values of (3. Increasing 
(3, as we have remarked above, means thickening the P tube, and we would therefore expect 
a heavier black hole to be required to expel this flux. This is indeed what we observe and, 
naturally, the curve for P is more affected by this than that for X. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

e e 

FIG. 4. Functions X(9) (left) and P{9) on the horizon for /3 = N = 1, a vacuum guess and 
M = 10, 2.5, 2, 1.9 and 1.8865. 




X P 
FIG. 5. Contours of X and P for M = Q = 1.8, N = (3 = l,e = 1(T 3 and N r = N e = 100. 

The same shift of the curves towards higher values of M happens when one increases N, 
as anticipated from the fact that this also thickens the string. Figure [7] shows the evolution 
of the critical mass M c (defined to be that at which the horizon cannot support a penetrating 
solution any more) in function of (3 (for N — 1) and of N (for (3 — 1). 

V. STRING ENDING ON A BLACK HOLE 

Next, we turn our attention to the case of a string ending on the black hole. This is 
an important configuration to consider, since it is the main 'phenomenological' input to the 
instantons mediating defect decay . Originally, [|12|] supposed that such a configuration 
may not be able to exist, however, the thin string arguments indicate that at least for large 
M, such a configuration is possible. What we will show is that while it is always possible 
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0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 

1/M 1/M 

FIG. 6. Values of X(tt/2) (left) and P(it/2) as functions of the inverse energy of the extremal 

black hole. This is done for N = 1 and the following values of 0: 500, 100, 20, 10, 5, 2, 1 and 1/2. 
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(3 N 
FIG. 7. Values of the 'critical mass' M c at which the only solution on the horizon corresponds 

to the string being expelled, in function of (3 (for N = 1) and N (for (3=1). 
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for a vortex to end on a black hole, for small M there is also a phenomenon analogous to 
flux expulsion: the X field is forced to sit in its unbroken phase (X = 0) on the horizon, 
and the P field takes the form of a monopole potential. 

To see this analytically, consider the horizon equations (0). These will have the bound- 
ary conditions X = 0, P = 1 at 9 = 0, and X = X m , P = at 9 = tt. The equation of 



motion for P, (|14b|) , can be integrated to show that 



A(l 



COS I 



<2P <(1 



COS I 



(28) 



where A = 1 — M 2 X^ n / (3. This shows that as M 2 //3 — > 0, P approaches its monopole form, 
P mon = |(1 + cos9). From now on, we will assume that M 2 / [3, M 2 /N 2 <C 1 and explore the 
possible solutions for X. First of all, note that X = 0, P = P mon is always a solution to the 
horizon equations, and it is this that we will call the expulsion solution. Now suppose that 
a 'piercing' solution exists, then ( |14a| ) implies that X has a local maximum at tt, which we 



will denote X m . Defining 9 by <9e(sin 9X^) = 0, we see that at 9q, Pq 



2N- 



sm 2 9 (l-X 2 ). 



Since M/N is assumed small, 9q will obviously be close to tt, and using the bounds on P, it 
is easy to see that 



2 1 n\ -a 4M 

— (TT — 9n I < S1I10 O < — F= 



(29) 



and so P will be extremely small. Integrating (|14a|) between and tt then gives 



X 



M 2 P 2 
sm9(l-X 2 )-- 



2N 2 



sin 9 



e 



X 



M 2 



sin 9 2N 2 



sin 0(1 -X 2 



(30) 



The LHS of this equation can be readily be bounded above using 



X 



M 2 
2N 2 



sin 9(1 -X 2 



P 2 
sin 9 



< 



4M 4 X„ 
N 4 X 2 



(31) 



The RHS is a little more tricky to bound below, but noting that X t ee is positive on 
[tt/2,9 ], yet negative (X^g = —M 2 X m {l — X^J/4) at tt, we can bound X t g on [tt/2,tt] by 



Xa < X m iV 2 sin 3 #n 



(32) 



and hence we can bound X below on (tt/2,tt) by 



X(9) > X v 



1 - 



AttM 3 
NX 3 



(33) 



therefore 



So 



X 



P 2 



M 2 



sin 9 2N 2 



sm9(l-X 2 



> I X 

/tt/2 



> x„ 



A^ 
16 



M 2 sin 9 



sin 9 2N 2 

4vrM 3 l / 8M 2 

Iva 1 " r 



(34) 
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Comparing these bounds on the RHS and LHS of ( |3~0|) we see that an expulsion solution is 
the only possible solution for 



A 4 A 4 
~6T 



M < 




8M 2 " 



(35) 



For N = (3 = 1 this gives M < 0.3. As before, this is a very weak bound, however, the 
important thing is that it shows that there is a lower bound on the values of M for which 
the X-field can vary on the horizon. For M sufficiently small, the Higgs field is forced to lie 
in its symmetric phase on the horizon, and we get an expulsion solution. 

On the horizon, the single string case differs from the one we have consider previously 
by the boundary conditions only. At 9 = 0, we must clearly have a string, but at 9 = n 
nothing forces the fields to assume a vacuum configuration. In fact, we have found that the 
only smooth solutions were such that X had a vanishing ^-derivative at the South pole. As 
figure H shows, the value of X m then depends on the black hole's mass, which means that 
if we wish to integrate the equations on the whole grid, we also have to update the 9 = rr 
boundary To find equations of motion on this line, we assumed that Pj sin 9 — > and that 
Xfi — > 0. The resulting equations, however, were particularly unstable against numerical 
errors. We were finally able to tackle this problem by artificially coupling the horizon to 
the rest of the grid (assuming continuity) for some hundred iterations, and by updating the 
fields at 9 = 7i assuming X t $ = Pg = there. Examples of a thin string solution piercing 
the horizon and of a thicker string being expelled are shown on figures |9] and 10|. Figure |8| 
shows the variation of the fields on the horizon as we thicken the string, illustrating that for 
M = 1 the Higgs field is already expelled from the horizon. 





1.5 2 2.5 3 0.5 1 1.5 

e e 

FIG. 8. Fields X{9) (left) and P(9) on the horizon for the case of a single string ending on the 
black hole. The string has a Higgs width of 1, and we plot the profiles for (3 = N = 1 and the 
following values of M: 9, 2, 1.3, 1.1, 1.03, 1. 



VI. GRAVITATIONAL BACKRE ACTION 



We begin this section with a lightning review of a self-gravitating cosmic string. For 
convenience, we will take N = 1 in this section. 
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X 

FIG. 10. Contours of X and P for a single 
[3 = N = 1, N r = N g = 100, r m = 10. 
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P 

ending on the black hole, and M = Q = 1/2, 
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We can readily extend the vortex to a self-gravitating system by using Thome's cylin- 
drically symmetric coordinate system [21 

ds 2 = e 2 ^dt 2 - e 2 ^ (dz 2 + dR 2 )- -a'e^d^ 2 (36) 



where 7, if), a are all functions of R only. 

In these coordinates, the energy momentum tensor becomes 



T» = £ 



„2</> y2 p2 p/2 
-2( 7 -V>) X /2 + e A_ + /3e 2( 7 -2V)£_ + I (X 2 _ 1)2 



-V, 



R 



e 2ip P2 p'2 
-2( 7 -V>) X /2 + e ^ - ^(7-2^ + I (X 2 _ 1)2 



„2V> \-2 p2 

-V 6 = e- 2 ^X' 2 - - pe 2 ^^— + \{X 2 - l) 2 



a" 



or 

p/2 

,2 ^ 4* 



a" 



T* = -V, = T n c 



To zeroth order in e = 87fGr] 2 introduced in section [IT] 
a = R , ^ = 7 = 0, X = X 
and conservation of the energy-momentum tensor gives 

(Rv oR y = v o4> . 

To first order in e the string metric is given by |22 



a 



1-e [ R R(£ -P oR )dR R + e [* R 2 (£ - V oR )dR 
Jo Jo 



7 = 2^ = e [ R RV oR dR. 
Jo 



(37a) 

(37b) 

(37c) 
(37d) 

(38) 

(39) 

(40a) 
(40b) 



Where the subscript zero indicates evaluation in the flat space limit. Since the string func 
tions X and P rapidly fall off to their vacuum values outside the core, the integrals in 
rapidly converge to their asymptotic, constant, values. Writing 



roo roo re 

A= R(S -V oR )dR, B= R 2 (£ -V oR )dR and C = 
Jo Jo Jo 

then the asymptotic form of the metric is 

ds 2 = e tC [dt 2 - dR 2 - dz 2 } - R 2 (l - eA) 2 e- eC d(j) 2 
= dt 2 - dR 2 - dz 2 - R 2 {1 - e(A + C)) 2 d<j) 2 , 

where t = e eC ^ 2 t etc. This is a conical metric with deficit angle 

A = 2ite(A + C) = 2vre J R£ dR = 8nGfi, 



RV, 



o-R' 



(41) 



(42) 



(43) 



where \i is the energy per unit length of the string. Notice that the deficit angle is inde- 
pendent of the radial stresses, but that when the radial stresses do not vanish there is a 
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red/blue-shift of time between infinity and the core of the string. The only case in which 
these stresses do vanish is when (3 = 1. 

Now we are ready to consider the gravitational effect of the string superimposed on the 
black hole. For this we need to consider a general static axially symmetric solution to the 
Einstein-Maxwell-abelian Higgs equations derived from the action S1 + S2, where Si is given 
by (0), and 5*2 is the Einstein-Maxwell action 



S2 = ~Y^q J d x ^9 [R + , (44) 

where for clarity in what follows, we have rescaled the electromagnetic field by \fG. We 
may pick coordinates so that the metric takes the form 

ds 2 = e 2 ^dt 2 - e 2 ^\dz 2 + dp 2 ) - a 2 e-^d(P 2 , (45) 

where V>7> a are independent of t, (p. [Note the deliberate similarities with the Thorne 
metric (0).] We then apply an iterative procedure to solving the field equations in the thin 
string limit, starting with the Reissner-Nordstr0m background solution, the Nielsen-Olesen 
forms of X and P, and expanding the equations of motion in terms of e = 8irGr) 2 as before. 

The usual Reissner-Nordstr0m metric (0) is of course written in 'spherical' coordinates, 
whereas, in such an iterative process we require it in axisymmetric coordinates. For future 
reference the coordinate transformation is 



p= ^r 2 - 2Mr + Q 2 sin 9, z = (r - M) cos 9 (46) 

and the metric and (rescaled) electromagnetic gauge potential (in a suitable gauge) are given 
by 

do 2 _ (fl + +fl_) 2 -4A 2 ,, 2 (R++R_+2Mr ( , 2 , 2 x 2 {R + +R_+2Mf , , 2 (AJ) 

Qb - {R++R.+2MY Ql 4R+R- \ a P ^ aZ ) P (R ++J R_)2-4A 2 ' V 1 ' ) 

, R ^l z R d v (j) magnetic 

A o»= S + 20 a . „i„,-„ (48) 



where 



i u±+R-+2M) electric, 



R± = (z± A) 2 + p 2 (49) 



and A 2 = M 2 — Q 2 . 

Returning to the general system, the relevant equations equations of motion are: 

= d p (ocF p v ) + d z (aF z n (50a) 

a zz + a pp = -V=9(V + V) ( 50b ) 

{a^, z ), z + {at P ), P = |V=#cr ° - t; - t; - r/) (so c ) 

7,pp + %zz = ~f P - VS - e 2 ^Tl (50d) 
where the energy momentum tensor is given by 
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T b a = El + eT b a . (51) 
The first term is the electromagnetic contribution to the stress-energy which is given by 



E$ = -F uX F>* + {F 2 5» 
and the last term is the contribution from the string, which has the explicit form 



X 2 P 2 

T ° = V(X) + 4 



ere 



-2i/> 



V(X) 



X 2p2 



+ 



TP 4- T z z = 2V(X) + 



2X 2p2 



(P 2 P + Pi) 
P~l a 2 e -2^ 

{Pi + PI 



HK + X l 
+ (x 2 p + x 2 z ) 



3 -2( 7 -V) 



-2( 7 -V) 



ere 



-2ip ' 



(52) 

(53a) 
(53b) 
(53c) 



Note that the electromagnetic stress energy always satisfies E p p + E z z = 0. 

As in we now write a = «o + e«i etc., and solve the Einstein-Maxwell and the string 
equations iteratively. The main difference to |J is that we now have the electromagnetic 
gauge potential present, which appears at 0(e°). This means that to 0(e) the geometry is not 
only affected by the string, but also by the backreaction of the string on the electromagnetic 
field. 

To zeroth order, we have the background solutions fl47f), (|48|), X = X (R) and P = Pq{R) 
where R = rsin# = pe~^° as before. Using the coordinate transformation (^6|), one finds 
that 



r2 z + R 2 P 



R+R- 



' 2M . Q 2 . " 

1 sin d H — - sm d 

ry ry £i 



3 2(7o-^o) 



(54) 



in the core of the string, where sin (9 = 0(M _1 ). Therefore, in and near the core of the 
string, the relevant combinations of the zeroth order energy momentum tensor of the string 
are 



7j = 5 + 0(M- 2 ) 
T*=-V^ + 0(M- 2 ) 



rpp _|_ rpz 

op ' oz 



(55a) 
(55b) 
(55c) 



As in || we will assume (and show subsequently that it is consistent to do so) that the 
perturbed solutions take the form 



a 1 = pa(R), if> 1 = il) 1 (R), ji = 7i(i?), 



.4 



f(R)A 



on ■ 



(56) 



Computing the necessary derivative of R gives from (|50b|) the following equation for a(R): 



r 2 (2M g 2 



a"(R) + 



2 Q 2 R (3M Q 2 



a'(R) 



2?"^ ^2 I ry ry^ 

-1 (W)' + 0(M- 2 ) = -{£ - V R ) + 0{M- 



(57) 
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which is consistently solved, as in the Schwarzschild case, by 

a(R) =- J R[S - V 0R }dR + - J R 2 [£ - V 0R )dR ~ -A + ^ as R -> oo. (58) 

Now look at the Maxwell equation for an electric potential (the magnetic potential can 
be obtained by a duality transformation). Substituting in the assumed form for the func- 
tions fl56|) , we obtain for f(R): 



f"(R) , /'(it!) p 2 



R 



+ 



R 2 



m 2 



(a'-2^) = 0(M- 2 ), 



(59) 



which implies that / = fo, a constant. Turning to the if) equation, and inputting in this 
form of /, we find 



1 + 



i? 2 /Q 2 2ikf 



+ 



R 



R 2 Q 2 



1 - 



+ 



Q 2 



a'R (Q 2 M 



+ 



\{V R + V< t) ) + ^ I [2U-2if) 1 



(60) 



This is solved by 



fh = \f RVr ~ | as i? -> oo, / = C/2 



(61) 



the latter value of fo being set by consistency of the if) equation outside the core. It is then 
straightforward to check that 71 = 2if)\. The magnetic correction is obtained either directly, 
or via duality, to be 



f M {R)=a{R)-2ip l {R) + C/2. 



(62) 



As in H, the corrections are almost the same as for the self-gravitating cosmic string. 
After transforming back to Schwarzschild coordinates, the metric outside the string core 
becomes 



ds 2 



cC 



(» 



2M 



+ ^)dt 



dr 2 



r 2 d6 2 



1 2M I Q 2 



r 2 sin 2 0(1 - eA) 2 e 



(63) 



where we have neglected the B term from a(R) since it yields a correction 0(Gfi) x 0(E x ). 
The gauge potentials are 



f Q(l - cos0)[l - e(A + C/2)] d u </> magnetic, 
{-£[! + eC/2]d„t electric. 



(64) 



If we now rescale the metric so that t = e eC ^ 2 t, etc, (and, accordingly, rescale the parameters 
M and Q) we find 



ds * = f x _ 2M + C£\ d p - f 2 d9 2 - f 2 sin 2 fl(l - eAfe 

\ r r2 J y — _|_ 01 



2 — 2eC j j 2 



(65) 
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Again, we find a deficit angle A = 2ire(A + C) = 8irGfi. Besides, the gravitational mass of 
the black hole, M g , which is given by the coefficient of 2f _1 in g%, has been shifted in the 
presence of the string to M g — M — e eC l 2 M. 

On the other hand, the internal energy of the black hole (its ADM mass, appropriately 
generalized to asymptotically locally flat spaces, see e.g., |23|1 ) is now 



Mi = M{1 - eA)e~ eC = M{1 - eA)e" eC/2 . (66) 

Not only is the mass of the black hole corrected, but the physical charge of the black hole, 
defined as 

q»= {Its r:r ^ 

becomes 

Q ph = Q[l-e(A + C/2)]. (68) 

Notice that M g /Q ph > Mi/Q ph = M/Q. 

We can now write the first-order corrected solution in terms of the physical parameters 
Mi, Q ph , /J,, as 

ds 2 =(l- e 4G ^ + e 8G ^)dt 2 — pg- - r 2 d6 2 - e~ 4G »r 2 sin 2 9dtf, 

f r 2 

(69) 



and 



J Q p h(l - cos8) d v <p magnetic, 
v ~ \ -e^Qf dj electric. [W) 

The corrected inner and outer horizons exterior to the core are therefore at 



r ± = e^(Mi±yjM?-Q 2 ph ), (71) 

since we are assuming that the string is thin with respect to the back hole, these expressions 
will hold except at the poles of the horizon. Notice that the condition that the black hole is 
extremal (i.e., its horizon is degenerate), is that Mi = Q p h- Finally, we can use the methods 
of to find the entropy using the Bekenstein-Hawking formula 

S = 4± = IQne^f 2 (72) 

Consider now an extremal black hole-cosmic string merger. If we keep the internal energy 
and the charge of the black hole fixed (i.e., impose microcanonical boundary conditions), 
then we find, as in 0, that the change in the gravitational mass, 5M g = 4G[iM g equals 
the energy of the length of string swallowed by the black hole. Since the entropy increases, 
we see that merging is thermodynamically favoured. And more remarkably, we see that, at 
least to first order, the black hole remains extremal after the merger. 
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VII. DISCUSSION 



Our analysis in this paper appears to settle the question of whether or not a vortex can 
penetrate an extremal black hole. We have provided analytical proofs that vortices of size 
smaller than a certain fraction (of order unity) of the radius of the black hole will definitely 
pierce the horizon, whereas vortices thicker than a certain lower bound will instead wrap the 
black hole. The numerical analysis confirms this, and, for N = (3 = 1, places the transition 
at M = 1.8865, or Gm = 1.8865/^. For a single string ending on the black hole the 
Higgs field presents a similar behavior, i.e., it vanishes on the horizon only if the black hole 
is small enough. In that regime, the magnetic field, instead of being expelled, takes the 
form of a monopole field. Thus we see that single strings are always allowed to end on black 
holes, which solves one of the paradoxes that the results of |T2j seemed to pose. Finally, 
we have computed the backreaction effect of a thin vortex on the geometry. This results in 
the expected conical geometry, but we have been able to check as well that the black hole 
remains extremal after including the corrections to the mass and the charge. 



Given that in [|T^] the expulsion of the (unbroken) magnetic field was related to a sort of 
'superconducting' behavior of the extremal black hole, one would be tempted to interpret the 
penetrating solutions as exhibiting the well-known breakdown of the superconducting state 



for strong enough magnetic fields. However, this does not seem to be the case here. In 
exact solutions (which account fully for the backreaction of the gauge field) are presented 
for extremal black holes in magnetic fields, and the expulsion persists no matter how strong 
the magnetic fields are taken to be. 

As a matter of fact, we can argue that, far from having anything to do with the strength 
of the magnetic field, it is instead the presence of a mass for the gauge vector which spoils 
the expulsion from the extremal horizon. In order to illustrate this point, consider a massive 
vector (Proca) field, with an explicit mass fi. On the extremal horizon the field equation 
becomes 

sm6d e f-^Yl -/i 2 P = 0. (73) 
\ sin 6 J 

This equation does not admit a constant solution (apart from the trivial P = 0). Therefore 
Ggcf, oc dgPcf, 7^ and we find a (locally) non-vanishing flux of the field across the horizon. 

As we have seen, things are subtler when the mass originates from spontaneous break- 
ing of the gauge symmetry. The system has both massive and massless phases, and both 
expulsion and penetration can be found, depending essentially on the relative values of the 
vortex and black hole radii. But the argument above shows that the expulsion of the field 
can take place only if the symmetry is exactly restored on the horizon. It is therefore quite 
remarkable that, in certain regimes, the geometry of the extremal horizon can locally enforce 
the exact restoration of the symmetry. 

The transition from penetration to expulsion can be viewed as a phase transition on the 
horizon of the black hole. In particular, Fig. ||] is reminiscent of the behavior of, say, the 
magnetization of a ferromagnet as a function of temperature, or more generally, the order 
parameter of a system undergoing a second order phase transition. The order parameter 
is in this case the value of the Higgs field on the horizon, X, and instead of a function of 
the temperature, the phase transition takes place when we vary the (inverse) size of the 
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horizon ~ M _1 . Beyond a certain critical value, M _1 > M" 1 , the symmetry is restored 
throughout the horizon. Notice that the transition takes place when the energy scale set 
by the black hole, M, is similar to the Higgs energy \f\r]. In a sense, this would be the 
natural expectation, since r\ sets the energy scale for symmetry restoration. But it should 
be stressed that this expectation is only realized for extremal black holes: the symmetry is 
never restored on non-extremal horizons]]. It might be interesting to pursue this analogy 
further, and study, e.g., critical exponents near the transition point, such as X ~ \M — M c \@ 
(and see, e.g., how j3 varies, or not, at different points on the horizon). 

Another interesting question that can be explored in more generality is the interaction 
between the straight cosmic string and the extremal black hole. The set up we have been 
considering so far places the vortex in perfect alignment with a black hole axis. If the 
black hole and the vortex are displaced relative to each other, the symmetry of the system 
decreases and complication increases greatly. We can, however, analyze in some detail the 
interactions between the extremal black hole and the cosmic string when the string is thin 
enough to allow us to effectively integrate out the details of the core structure. We can 
proceed in several levels of approximation. A very crude approximation would be taking the 
black hole as a test particle in the background of a self-gravitating cosmic string, i.e., in a 
flat spacetime with a conical defect. As is well known, since the spacetime is locally flat the 
test particle does not experience any force. 

We can improve on this by accounting for the gravitational field of the black hole. In 
the Newtonian approximation, we would be solving the Poisson equation for the Newtonian 
potential in a conical spacetime. The effect of the conical defect on a massive neutral particle 
can be understood by viewing the particle as subject to a force coming from the 'images' 
produced by the conical defect. This results in the neutral particle being attracted towards 
the string [TBj (see also pi|). In contrast, when applied to a charged 'extremal particle,' 



such that m = q, this argument would yield a vanishing force, since the gravitational and 
electrostatic forces between the particle and its images cancel out. 

Actually, for extremal black holes this result holds not only in the Newtonian approx- 



imation, but also exactly in the full Einstein-Maxwell theory |25]. To see this, recall first 



that the solution for an (electric) extremal black hole can be written as 

ds 2 = H~ 2 dt 2 - H 2 (dx 2 + dy 2 + dz 2 ), (74) 
A M = (H- 1 - 1)8^. 

Written in this fashion, the Einstein-Maxwell equations only require that H is a harmonic 
function in the flat (x,y,z) space, V 2 xy z H = 0. The extremal black hole is recovered by 
setting 

H = l+ Q (75) 

yj(x- x ) 2 + {y- y ) 2 + {z- z ) 2 



Notice as well that it would be incorrect to think that the restoration of the symmetry comes 
about as an effect of the thermal properties of black holes. For one thing, an extremal black hole 
has zero temperature. Moreover, the thermality is only seen when accounting for quantum effects, 
whereas here we work at the classical level. 
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and the horizon is at (x,y,z) = (xq, yo, ^o)-0 If we want to include the cosmic string, then 
we just have to solve the Laplace equation for H, this time in a space with a straight conical 
line. The relevant solution has been given in many places, see e.g., in the context of cosmic 



strings, |1| . In cylindrical coordinates (p, z, 0) centered on the string, with conical deficit 



such that < (f) < 2tt/p, if we put the black hole at p = p , <p = and z = 0, then 

Q [™ du psinhpw 

Hp{z, p, <p; po) = 1 + — /== / . == — r -, 76 

7i^/2ppo Ju vcoshw — coshwo coshpw — cos pep 

where Uq is defined by 

eosh„ ^ 2+ / + "° . (77) 

This solution is nonsingular, away from the conical line and the singularity of the black hole. 
Thus, there are no forces between the extremal black hole and the cosmic string. 
At distances much larger than p , the harmonic function becomes 

H *- l + WT?- (78) 



Since p ~ l + 4Gp>, we see that we reproduce the gravitational and electric potentials in (p9|), 
([70|), with the parameter Q equaling both the internal energy Mi and the physical charge 
Qph w e had introduced earlier. The result, however, is independent of whether the string 
and the black hole are merged or not. Whenever the black hole is in the presence of a 
cosmic string, the lines of force are 'focused,' resulting in an increase of the gravitational 
and electric potentials at long distances. 

The same conclusion regarding the absence of a force between the two objects can be 
reached again from another perspective. This time, neglect the gravitational backreaction 
of the string, and consider a Nambu-Goto string in the background of the extremal black 
hole, with metric as in (|74]) . The static interaction potential V(x,y,z) experienced by a 
string can be read off from the string action, Ing ~ / drV(x,y,z). If we take a straight 
string along, say, the z axis, this is, T = r, Z = a, X, Y = constant, then the action is 

Ing = -p< J dr da^J- det g liy d a X^d p X u 

= -pL J dr , (79) 

(L is the length of the string) i.e., the static force vanishes. With a bit more work it is easy 
to see that if the string is given a velocity transverse to its axis, then its motion is slowed 
down as it approaches the black hole. 

Clearly, in all these arguments we have been neglecting the effect of the black hole on 
the string core. But nevertheless we seem to find that, at least if they are well separated, a 
straight string and an extremal black hole will hardly feel each other's presence. 



2 The 'Sell war zschild coordinates' used in (12) correspond to setting r — M = [{x — xq) 2 + (y 
y ) 2 + (z - ^o) 2 ] 1/2 . 
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These results appear to imply that the binding energy between the infinitely thin string 
and the extremal black hole is zero. It appears difficult to make a clear comparison between 
the energy of a finite-radius vortex before and after the merger. On the one hand, when 
the vortex penetrates the horizon, a hole is cut out from space, where a part of the string is 
missing. As we have seen in the previous section, the energy of the missing link equals the 
change in the gravitational mass of the black hole. On the other hand, since the geometry 
along the string changes in the process, one should be careful about how to define the total 
energy change. In particular, since the string has an infinite length (and hence, infinite 
energy), one should specify a regularization, and choose how to fix a large but finite length 
of the string before and after the merger. 

We conclude by mentioning that it should be interesting to include fermion fields and 
study the supersymmetry of the extreme black hole-vortex configuration. As is well known, 
in the absence of the cosmic string, the Reissner-Nordstr0m black hole can be embedded in 
(N = 2) supergravity, and the extremality condition M = Q appears then as the BPS con- 
dition for the existence of unbroken supersymmetry generators ||26|| . Similarly, the Nielsen- 
Olesen vortex with (3 = 1 admits a supersymmetric embedding in N=2 supergravity in three 
dimensions |27| , and the solution preserves half of the Killing spinors of the flat vacuum. A 
natural question to ask is whether the merger configuration will be supersymmetric as well. 
Although we have not analyzed this point in any detail, our analysis indicates that when the 
vortices are infinitely thin, the system exhibits some of the features characteristic of BPS 
configurations, including the equality of the ADM mass and the physical charge, Mj = Q p h 
(which was actually obtained for finite, if small, width strings), as well as the absence of 
forces and vanishing binding energy between the two objects^. Even when the vortices have 
a finite width, we suspect that it is also possible to maintain another feature of BPS sys- 
tems, namely the reduction of the equations of motion to a first order system. Although at 
first sight this seems unlikely given the lack of symmetry, it was shown in Anderson et. al. 
P8fl that for a general worldsheet embedding, such a reduction to a first order system does 
occur. In this case we have a curved geometry, nonetheless, preliminary indications are that 
a generalization of this method will work. However, the drastic change in the behavior that 
takes place as the vortex grows thicker (from penetration to expulsion) would make it very 
surprising that supersymmetry be present in general. 
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